Boundary hopping and the mobility edge in the Anderson model in three dimensions 
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It is shown, using high-precision numerical simulations, that the mobility edge of the 3d Anderson 
model depends on the boundary hopping term t in the infinite size limit. The critical exponent v is 
independent of it. The renormalized localization length at the critical point is also found to depend 
on t but not on the distribution of on-site energies for box and Lorentzian distributions. Implications 
of results for the description of the transition in terms of a local order-parameter are discussed. 



Since the discovery of the Anderson localization [lj], 
one of the central problems of the physics of disordered 
materials was to understand the metal-insulator tran- 
sition induced by the localization of electronic states. 
The scaling theory of localization proposed that the 
disorder-averaged dimensionless conductance (g) is the 
only relevant scaling parameter governing the transi- 
tion, and that there is only one fixed point under the 
renormalization-group transformation (RG) in Id and 2d 
corresponding to the insulating phase, with 2 being the 
lower critical dimension of the transition, while in 3d 
there are three such points and a continuous localization- 
delocalization transition characterized by a critical expo- 
nent v. 

The discovery of the universal conductance fluctua- 
tions Q showed that g is not a self-averaging quantity, 
and therefore scaling of its whole distribution has to be 
studied. Nonetheless, there is a compelling evidence that 
the scaling properties of the distribution itself in the crit- 
ical region of the transition are still governed by a single- 
parameter [3]. 

The transition has a field-theoretical description within 
the super-symmetric nonlinear er-model [5] (NLS), but 
technical difficulties with the calculation of v from the 
2 + e expansion do not yet permit accurate determina- 
tion of v, despite significant theoretical progress [f|. v is 
currently most accurately determined numerically, espe- 
cially via the transfer-matrix method (TMM) 0,0, lUEj], 
also used in this work and described below. TMM gives 
the phase diagram of the transition, by calculating the 
critical disorder strength, W c , at which a critical state 
appears at the band center {E = 0). For W < W c there 
is a band of extended states within \E\ < E C (W). The 
mobility edges W c and E C (W) depend on the distribution 
of impurity energies [IH . 

The subtlety of the Anderson transition is that natu- 
ral choices for an order parameter, such as the density of 
states or the Q fields, that are conjugated to the symme- 
try breaking field in an effective field-theoretical descrip- 
tion of the transition within NLS, behave analytically in 
the vicinity of the transition 0, [l3]. Nevertheless, a de- 
scription of the transition in terms of an order parameter 
function (OPF) has been proposed [l3[ and shown to be 
related in a deep way to the distribution of the local den- 
sity of states [14l | . 

Another surprising feature of the transition is the in- 



fluence of boundary conditions (b.c) in the critical region 
on several important quantities: critical distribution of 
g [HI, critical level spacing E3, EH, critical renormal- 
ized localization length A c [19] (that is also studied in 
this work) as well as on the /3-function of (g) in diffusive 
approximation [la ]. These are generally understood in 
terms of the dependence of the quantized diffusion modes 
on the choice of b.c [HI, the sample geometry [2(| and 



topology [2l|. The case of general boundary conditions 
was recently considered in Ref. [221 ]. 

To further clarify the role of b.c and surface effects 
on the properties of the localization-delocalization tran- 
sition, a model is proposed here that corresponds to the 
Anderson model with boundary hopping terms taken to 
be continuous variables. In this model one expects diffu- 
sion modes to be quantized in a quantitatively different 
way than in the models with periodic or open b.c. The 
principal finding of this Letter is that the mobility edge 
of the Anderson model in 3d depends on the boundary 
hopping term. 

The studied Hamiltonian is: 
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where i,j denote sites of the simple cubic lattice of 
L x x L y x L z atoms; |i) is the eigenstate of the electron 
localized at the i-th site with impurity energy ej, which 
is an i.i.d random variable with the probability distribu- 
tion P(e,W), here taken to be box distribution (BD), 
i.e. ti are uniformly distributed in [— W/2, W/2), and 
Lorentzian distribution (LD), P(e, W) = W/tt(W + e 2 ). 
The hopping integral in the bulk, tb, is set to 1. By 
taking t XtVlZ to be or 1 one obtains, as special cases, 
models with open, periodic and mixed boundary condi- 
tions studied in Ref. 18, lj| . 

This work studies the geometry of quasi-ld slabs L = 
L x = Ly and L z = M 3> L, (the density of boundary 
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hoppings is thus L~ v ) : and the localization properties 
of the 3d samples are deduced from the scaling proper- 
ties of the smallest Lyapunov exponent when L —> oo. 
To this end the eigenproblem H^f = E 1 ^ is rewritten 
in terms of the transfer matrices, \Pm = Tm^o, Tm — 

Tlf=i ( ^ l i ^ 0^1' wnere Hi'^i are ' respectively, 
projections of H on the i-th slice of the slab. The 
product converges towards the limiting symplectic ma- 
trix when M — > oo, with corresponding 2M 2 Lyapunov 
exponents 7$ being eigenvalues of limM^oo 53? m ^jv/'^A/j 
which are calculated to a given relative accuracy by mul- 
tiplication of a sufficient number of % combined with the 
Gramm-Schmidt reorthonormalization [1, [I, [23[ . Fermi 
energy is set E — 0, since the critical states are expected 
to appear first at the band center when W is decreased. 

The renormalized inverse of the algebraically smallest 
7i, A = (Lrnin^ | T« | ) 1 , defines the largest length scale in 
the system that is identified with the correlation l eng th, 
the scaling law of which near the critical point is [24j : 



A (to, L) = F(ip(w)L 1 ^,(p(w)L v ), w 



W-W c 



(1) 



where v > is the critical exponent, y < is the irrele- 
vant critical exponent due to the finite-size correction to 
scaling and W c is the critical disorder strength. Near the 
critical point (w = 0, L^ 1 = 0) the scaling function F, 
the relevant scaling field ip and the irrelevant scaling field 
(j> are expected to be analytic [24[ , and each of them can 
be expanded giving the following class of model functions 
for the fit: 



A(t«,L) = F (ip(w)L 1/l/ ) + 



(2) 



F n {ip{w)L 1/v ) 



]T <p(w) n L ny F n (i/;(w)L 1 / 1 '), (3) 
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(4) 



m=0 



, <j>(w) = fa™ 3 ; ( 5 ) 

3=1 j=0 



where coefficients of expansions are fitting parameters 
(except for Fqi — Fio — 1 because two parameters 
are not independent [l0|), that are determined from the 
least-x 2 nonlinear fit. We generalize the class of model 
functions from Ref. [1(3] by allowing each F n to be ex- 
panded up to fc n -th power, where k n are not necessarily 
all equal, and by allowing some of <f>j to be set to for 
j < ran, giving Tn' R < tur coefficients of the expansion 
of tp(w). The total number of fitting parameters is thus 
N P = En= (kn + 1) + m' R + mj - 1. 

Fq is the one-parameter scaling function, while Eq. ([3]) 
represents corrections to scaling due to the finite-size ef- 
fects described by y and (p. The corrections to scaling 
were studied in the context of the in teger quantum Hall 
effect [H, [H, H3] - They were shown fati^ to be impor- 
tant for the accurate estimation of v, while the formu- 
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TABLE I: Input parameters of the fit and values of obtained 
X 2 and Q. Nd is the number of points, 
fitting parameters. In all cases = ko 
Parameter tfe is set to for t = 0, 0.1. 



N p is the number of 
= 3, ni — 1, m/ = 0. 



lation of Ref. [l(| allows also accurate determination of 
W c . 

The case t = t x = t y is studied first. To determine the 
dependence on t, Lyapunov exponents were calculated 
and the scaling analysis done independently for several 
values of t. A relative accuracy of 0.05% was chosen 
in all cases studied, and correspondingly the slabs were 
up to 4 x 10 7 sites long. Values of N p fitting parame- 
ters {oii} = {y, y, W c , F mn , ip n , <p n } were determined us- 
ing the Levenberg-Marquard algorithm. 95% confidence 
intervals of are determined explicitly, by calculat- 
ing projections of the confidence region X 2 ({ a i}) < Ay 2 , 
where A\ 2 depends on N p and the confidence level [2a ] . 
This is a nontrivial calculation for which an efficient al- 
gorithm was developed. 

Table 1 summarizes parameters of the simulation, ob- 
tained values of \ 2 an d the "quality of fit" parameter Q. 
ft was suggested in Ref. [l(| that acceptable fits should 
have Q > 0.1, which the results presented here satisfy. 

The calculated values of A(W, L) for t = 0, 0.25 and 
1 are presented in Fig. 1, together with the lines cor- 
responding to the fits. The Fig. shows an increase of 
A(L) with L for smaller values of W and a decrease for 
larger W, which correspond, respectively, to the delocal- 
ized (metallic) and localized (insulating) behavior. The 
same values with the corrections to scaling, Eq. (J3]), sub- 
tracted are shown in Fig. 2, and the critical point is where 
all curves cross. Figure 3 shows the data collapse, where 
the upper branch corresponds to the metallic (w < 0) 
and the lower branch to the insulating (w > 0) phase. 

The obtained values of v, A c ,y, and W c are presented, 
respectively, in Fig. 4 and 5. The results in the case of 
open (t = 0) and periodic (t = 1) boundaries for BD are 
in excellent agreement with Ref. 13, 13] ■ They support 
the universality of v, and non-universal, i-dependent, A c , 
W c and y. A c (t) is, however, found to be independent of 
the distribution of e, which therefore supports the uni- 
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FIG. 1: (Color online) A(W, L) calculated for t = 0,0.25, 1, 
going from the bottom set of points to the top, with corre- 
sponding fits (lines). In each set the fastest changing A(W) 
corresponds to the largest L studied. Sets of points for inter- 
mediate values of t smoothly interpolate between the values 
in the Figure and are omitted for clarity. 



versality of the function A c (t). 

Most remarkably, the results suggest that W c (t) is not 
constant, i.e, the boundary hopping term t induces a shift 
of the critical point of the infinite system. To understand 
this better, we recall that in the classical systems with 
short-range interaction and the local order-parameter, 
one generally expects surface ordering with a different 
critical temperature T c from that in the bulk, but the 
T c of the bulk is not affected by the presence of bound- 
aries |24j . If these assumptions are generic for any theory 
with local order-parameter and short-range interactions, 
and having in mind that disorder plays the analogous role 
to temperature in classical phase transitions, one arrives 
to the conclusion that there is no such theory in the uni- 
versality class of the H studied here for generic values 
tx,y,z 7^ 0, 1. It should be noted that the non-constant 
W c (t) found here is not in contradiction with the OPF 
description of the transition, since OPF at a given point 
is defined as an integral involving values of Q-fields at all 
other points of the sample. 

In addition to H with constant t, systems where each 
individual boundary hopping term, ty, is an i.i.d random 
variable were studied for three distributions: tij 6 [0, 1), 
Uj G [—1, 1), and when Uj takes the values or 1 with the 
probability 1/2. In all three cases BD is used for a. The 
results for v, W C ,A C and y in all three cases give values 
close to those obtained for t — 0.5, thus suggesting that 
randomness in boundary hopping terms can be described 
by H with an effective constant t. 



FIG. 2: (Color online) Same as in Fig. 1 after the subtraction 
of the correction to scaling due to the irrelevant field. 
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FIG. 3: Data collapse for t = 0, 0.25, 0.5, 0.7, and 1 (from the 
bottom two curves to the top). 

Furthermore, the t x ^ t y case was studied as well for 
t x = 0.7, t y = 0, when A c w 0.47(2) was obtained for 
both BD and LD. Since they differ from A c (t x — t y — 
0.7) — 0.543||, it seems reasonable to expect that there 
is also a universal A c (t x ,t y ) function independent of the 
distribution of on-site energies. 

Finally, the t = —0.1 case with BD of e$ was stud- 
ied to see effects of the sign of t on critical properties, 
and results for v, W c , A c , and y are approx. the same as 
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FIG. 4: (Color online) Obtained values of the v, A c and y for 
box (filled circles) and Lorentzian (open circles) distributions 
of disorder. The latter are slightly displaced along the x-axis 
for clarity. Error bars are 95% confidence intervals. 



16.8 



16.7 



16.5 




4.35 ■ 



Lorentzian 



0.0 



0.5 
t 



1.0 



FIG. 5: (Color online) Dependence of the critical disorder 
strength W c (t) (mobility edge) for the box (left panel) and 
Lorentzian (right panel) distribution of on-site energies, sepa- 
rating the metallic (W < W c ) from the insulating (W > W c ) 
phase. The value W c (t = 1) in the Lorentzian case is from 
Ref. Error bars are 95% confidence intervals. 



those obtained for t = 0.1. This is not surprising since 
it was already found [l9[ that the change from periodic 
to antiperiodic b.c. does not affect critical properties. 
From these two results it seems reasonable to infer that 
Ac(t) w Ac(-t). 

In conclusion, the Anderson model with the bound- 
ary hopping term t is studied near the localization- 
delocalization transition in the presence of the time- 
reversal symmetry for two distributions of on-site ener- 
gies. After the careful estimation of statistical errors, the 
universality of v and A c (i) is shown, and found that W c 
depends of t, which is interpreted as the absence of the 
local order-parameter description of the transition. The 
surprising deviation of up to 1% of W c (i) remains to be 
explained, for example by a theory that would account 
for the occurrence of the maximum of W c at t « 0.25. 
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